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A mathematical  model  for  predicting  the  performance  of 
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discussed  in  detail  and  then  applied  to  a number  of  typical 
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PREFACE 


In  order  to  quickly  assess  the  effect  of  certain  types  of 
terrain  irregularities  on  the  perforaance  of  1LS  Glide  Slope 
antenna  systems,  a mathematical-electromagnetic  scattering  coaputer 
aodel  has  been  developed. 

This  work  was  perforaed  by  members  of  the  Modeling  and 
Analysis  Section  of  the  Electromagnetic  Technology  Division, 
Transportation  Systeas  Center  for  the  Category  1 and  2 Il-S  Section 
of  the  Terninal  Navigation  Branch  ol  the  Systeas  Research  and 
Developaent  Service,  Federal  Aviation  Ada ini strati on  (FAA) . 
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PART  I.  THEORY 

1.  INTRODUCTION 

A mathematical  model  has  been  developed  for  predicting  the 
performance  of  image-type  glide  slope  arrays  in  the  presence  of 
certain  types  of  terrain  irregularities.  The  basic  theory  is 
developed  in  the  following  section  and  then  numerical  results 
are  presented,  illustrating  the  performance  of  different  glide 
slope  array:  for  several  terrain  configurations. 
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2.  THEORY 


The  vertical  antenna  pattern  of  an  image-type  glide  slope 
array  is  determined  by  the  interference  between  the  direct  and 
ground- ref lected  radiation.  The  ground-reflected  radiation  is 
of  course  generated  by  the  charges  and  currents  induced  in  the 
ground  plane  by  the  incident  ndiation  from  the  glide  slope 
array.  For  simplicity,  we  will  assume  that  the  ground  is 
perfectly  conducting.  For  a perfect  conductor,  the  surface 
current  density  I is  given  by: 


1c  » ft  x ft. 


(1) 


where  ft  is  the  unit  normal  yector  pointing  out  of  the  ground 
and  ft  is  the  total  magnetic  field: 


H - Hi  ♦ Hs, 


[2) 


where  H|  is  the  incident  and  ft?  is  the  scattered  magnetic 
field. 


In  terms  of  the  surface  ■ urrent  density  t,  the  ground- 
reflected  or  scattered  magnetic  field  Hs  at  a receiver  above 
the  ground  is  given  by  the  following  surface  integral: 

Hs(Tj)  - J !(r)  x vGfrlcf)ds  (3) 

S 

The  position  vectors  r^  and  r denote,  respectively,  the  position 
vector  of  the  receiver  and  the  position  vector  of  a source  point 
on  the  ground  plane  S relative  to  some  arbitrarily  chosen  origin 
of  coordinates.  The  two-point  Green's  Function  G(r^,r)  is 
given  by: 


G(rltr) 


.iklVrl 

i ri *r| 


(4) 


where  k * ls/\  and  \ is  the  wavelength  of  the  incident  radiation. 
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We  will  adopt  the  physical  optics  model  for  the  surface 
current  density  Specifically,  we  will  assuse. that  on  those 
portions  of  the  ground  plane  not  directly  illuainated  by  the 
glide  slope  array  (shadowed  regions),  t is  identically  zero 
and  that  on  those  portions  of  the  ground  which  are  illuainated, 
t is  given  by  twice  the  tangential  component  of  the  incident 
aagnetic  field 

h 0 on  S.;  (S) 

t * 2 (A  x fi.)  on  St,  (6) 

where  S_  and  S+  denote,  respectively,  the  unilluoinated  and 
illuainated  portions  of  the  ground  plane  S.  The  physical 
optics  approxiaation  for  the  current  distribution  should  be 
reasonably  accurate  when  the  characteristic  dimensions  of  the 
terrain  irregularities  are  large  compared  with  the  wavelength 
X.  Substituting  equations  (5)  and  (6)  into  equation  (3),  we 
obtain  the  following  approximate  expression  for  the  ground  - 
reflected  magnetic  field  at  the  receiver  point  r^ : 

Hs(ri}  - ^ f lft(r)  x \(r))  x fc(?x,r)ds.  (7) 

S+ 

The  numerical  evaluation  of  the  integral  in  Equation  (7) 
for  a ground  plane  exhibiting  arbitrary  terrain  variations  is, 
in  general,  a prohibitively  time -consuming  operation.  Con- 
sequently, we  will  confine  our  attention  to  terrains  which 
vary  only  along  one  cocrdinate  axis.  Specifically,  we  will 
treat  ground  planes  which  exhibit  irregularities  whicn  do  not 
vary  in  the  direction  perpendicular  to  the  centerline  of  the 
runway.  The  touchdown  point  on  the  centerline  of  the  runway 
opposite  the  glide  slope  array  will  be  chosen  as  the  origin  of 
coordinates.  The  z-axis  is  chosen  to  be  the  vertical  axis 
passing  through  the  origin  0 and  pointing  out  of  the  ground 
while  the  x-  and  y-axes  lie  in  a horizontal  plane  whose  normal 
parallels  the  z-axis.  The  x-axis  points  along  the  centerline 
of  the  runway  while  the  y-axis  is  perpendicular  to  the  center- 
line. The  terrain  irregularities  are  assumed  to  vary  only  with  x. 


That  is  to  say,  the  equation  of  the  surface  of  the  ground  is  of 
the  fora  z - f(x)  and  is  independent  of  variations  in  y.  To 
put  it  another  way,  the  local  outward  noraal  to  the  ground,  fi, 
has  only  x and  z components. 

We  will  now  consider  the  problea  of  a horizontal  dipole 
nounted  above  such  a ground  plane  and  oriented  parallel  to 
the  y-axis,  the  invariant  axis  of  the  terrain  irregularities. 

The  geosetry  of  a typical  problea  is  illustrated  in  Figure  1. 

The  y-axis  points  into  the  page  and  the  entire  ground  plane  can 
be  generated  by  projecting  the  cross  section  depicted  in 
Figure  1 froa  -*  to  ♦«*  in  the  y-direction. 

To  calculate  the  scattered  field  Hs,  generated  by  a given 
terrain  contour,  the  ground  plane  is  first  divided  into  a nuaber 
of  connected  planar  sections  which  extend  froa  -»  to  ♦«*  in  the 
y-direction.  The  side  view  of  such  a section  is  represented  by 
the  darkened  segaent  of  the  contour  depicted  in  Figure  1.  The 

points  (xoi»:oi^  and  (x0,,202^  rePrescat  the  end  points  of  the 
section.  The  length  L of  a given  segaent 

(L  • V<x02-*01>2  * ^OZ'W2) 

is  so  chosen  that  the  unit  noraal  ft  is  essentially  constant 
for  z01<z<z0-,  and  xgjfxfXg,.  0nce  t*,e  contour  has  been  segnented  in- 
to these  planar  sections.  Equation  (7)  :an  be  applied  to  each  sec- 
tion and  the  total  scattered  field  can  be  obtained  by  suaaing  the 
contributions  froa  all  the  sections.  We  will  therefore  concentrate 
upon  obtaining  an  expression  for  the  contribution  to  the  total 
scattered  field  Hs  froa  a typical  planar  segaent  such  as  the 
one  depicted  in  Figure  1.  It  will  be  assuaed  that  only  those 
segaents  which  are  in  front  of  the  receiver  (-«<x<Xj3  contribute 
to  Hs-  The  terrain  behind  the  dipole  (-»<x<G)  is  assumed  to  be 
perfectly  flat.  The  exact  nature  of  terrain  behind  the  antenna 
is  not  really  critical  since,  in  practice,  the  dipoles  which 
make  up  glide  slope  arrays  arc  screened  to  eliminate  the  back 


course. 


Geometry  of  a Typical  One-Dimensional  Ground  Reflection  Problem 


Let  ft  denote  the  contribution  to  the  total  scattered  field 
s 

from  the  planar  section  shown  in  Figure  1: 

W " 2*  f IfixHjC?)]  x?G{r1(f)ds  (8) 

where  the  integral  is  carried  out  over  the  surface  of  the  segment 
(cf.  Eq.  (7)).  Let  §x>  £y,  and  6z  denote,  respectively,  the 
unit  vectors  in  the  x,  y,  and  z directions.  In  terms  of  these 
vectors,  the  unit  normal  n is  given  by: 

n * -sine  ex  + cose  e2,  (9) 

where  e is  the  angle  between  the  planar  section  and  the  x-axis 
(cf.  Fig.  1).  The  (x,z)  coordinates  of  the  transmitting  dipole 
and  the  receiver  are  denoted,  respectively,  by  (0,h)  and  (Xj,Zj). 
For  simplicity,  we  will  assume  that  the  receiver  is  on  the  center- 
line  of  the  runway  so  that  y^  * 0.  The  y-displacement  of  the 
dipole  will  be  denoted  by  y . Glide  slope  arrays  are  generally 
located  several  hundred  feet  from  the  centerline  of  the  runway. 

The  incident  magnetic  field  H^r)  produced  by  the  transmitting 
dipole  is  given  in  terms  of  the  volume  current  density  3 in  the 
dipole  by  the  following  formula: 

H.(r)  = Jy  x ?*G(r,r*)jdv\  (10) 

where  V*  operates  on  the  primed  coordinates.  Although  the 
actual  dipoles  which  make  up  glide  slope  arrays  are  half-wave- 
length bent  dipoles,  we  will  treat  the  case  of  a point  dipole 
with  the  following  volume  current  distribution: 

3(r')  » JQ  ey  6(x’)  6(y'-ya)  6(z'-h).  (11) 

This  simplifying  assumption  expedites  the  evaluation  of  the 
integrals  (10)  and  (8)  and  should  not  be  too  restrictive  since 
we  are  only  interested  in  the  vertical  radiation  pattern  formed 
by  the  direct  and  ground  reflected  radiation.  In  Equation  (11), 

JQ  is  some  constant,  the  6-functions  are  the  Dirac  delta 
functions,  and,  as  previously  noted,  y_  and  h denote,  respectively, 
the  displacement  of  the  dipole  from  the  centerline  of  the  runway 
and  the  height  of  the  antenna  above  ground. 

2-5 


”1 


Substituting  Equation  (11)  into  Equation  (10)  and  integrating, 
we  obtain  the  following  expression  for  the  incident  magnetic  field 
H^r): 


_ ^ kJ  ikD 

fi Ar)  - i-5-  £ 

1 4 it  D 


f(D)  [x  dz  - (z-h)gxJ, 


(12) 


where  r * x 8 ♦ y 8 ♦ z 8 and  D is  the  distance  from  the 

* y i 

antenna  to  the  field  point  r: 


(x2  ♦ (y-ya)z  ♦ (z-h)2]1/2. 


The  function  f(D)  is  given  by 

“D>  " B * E52- 


(13) 


(14) 


It  will  be  noted  from  the  form  of  f(D)  that  we  are  retaining 
the  near  field  as  well  as  the  far  field  contributions  to  ft.. 

Equations  (9)  and  (12)  can  now  be  substituted  into  Equation  (8), 
and  the  magnetic  field  Ks  at  the  receiver  point  = Xj  6^  + Zj  ez> 
due  to  the  single  planar  section  of  terrain,  evaluated.  Ke  will 
assume  that  the  receiver  is  in  the  far  field  of  the  terrain 
segment  so  that$G(rj,r)  can  be  written  as  follows: 

(r.-f)  eik'VfI 

7G(r.,r)  = -ik . f 1 51 

X , + I I “•*  -*•  ! 

I r1-r | 1 r j - r | 

Actually,  only  the  z-coraponent  of  Ks  need  be  calculated  since 
the  receiving  antenna  responds  primarily  to  this  component 
(equivalently,  the  horizontal  component  of  the  electric  field). 

The  z-component  of  Ks  will  be  denoted  by  hgz.  Substituting 
Equations  (9),  (12),  and  (15)  into  Equation  (8),  we  obtain  the 
following  expression  for  hgz  at  the  receiver: 


hsz(rl} 


kZJ0  f(xx-x)(x  sine  -(z-h)cose)f  (D)eik(1)+R)  Js> 

S,2  J 


DR 


(16) 


where  R * [(Xj-x)2  ♦ y2  ♦ (z^-z)2]1/2  js  jistancc  froo  a point 
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on  the  terrain  segment  to  the  receiver.  The  element  of  area  ds 
is  equal  to  dydn,  where  n is  a variable  ranging  from  0 to  L 
along  the  surface  of  the  terrain  segment  perpendicular  to  the 
y-axis.  The  x-  and  z-coordinates  of  a point  on  the  surface  of 
the  planar  section  can  be  expressed  in  terms  of  n as  follows 
(cf.  Fig.  1): 

x « xft1  ♦ r,  cose,  0<n<l. 

01  * • (17) 

z = Zqj  + n sine,  0<n5L 


Substituting  Equations  (17)  into  Equation  (16),  we  obtain 


tJ 

J0  -<*> 


- k2jo 

hsz(rl^  = 7T  (U0rh)cose  - x0l  sin£l  X 

(X,-Xfli-n  cose)  ik  fD+Rl 

1 01  f(D)  eliclu  KJdydn, 


OR 


r 


(18) 


where  D and  R are  expressed  as  functions  of  the  variables  y and 
n as  follows: 


D = [(x01<r,cose)2  ♦ (y-ya)2  ♦ (z01 -h^nsine) 2 J 1/2,  (19) 
R * ( (Xj-XQj-ncose)2  ♦ y2  ♦ (Zj-Zpj-nsine)2] 1/2.  (20) 


The  integration  with  respect  to  y(-»<y<+®)  can  be  performed 
approximately  using  the  method  of  stationary  phase  leaving  only 
an  integral  with  respect  to  n to  be  evaluated  numerically: 


h5Z(ri>  " 


k2J 


8 TT 


0 ((zQ1-h)cose  - xnisincjXi/,2ei  1,(21) 


‘01 


where 

L 

I » f e’ kF^n^G(n)dn , 

Jo 

and  F and  G are  defined  as  follows: 


(22) 


A - ( (xj-XQ.-’-cose)2  ♦ (Zj-Zgj-nsine)2]1^2. 

8 * ( (x01>ncosc)2  ♦ (z01-h>nsin€)2l1/,2t 

r y2  1 1/2 

C - 1 ♦ —2 — J , D * CB,  R - CA, 

L (A+B)  2J 

•■Hi -S’ 

F(n)  « C(A  ♦ B), 

f(D)  [x.-x^.-ncosc] 

G(n)  - i— ^ . (23) 

dr2s1^2 

The  contribution  to  the  total  scattered  field  Hs  from  the 
region  behind  the  antenna  (-®<x<0)  can  be  calculated  approximately 
using  asymptotic  methods.  As  noted  previously,  the  terrain  behind 
the  glide  slope  array  is  assumed  to  be  perfectly  flat  (e=0,  z*0 
for  -®<x<0).  From  Equation  (16),  the  contribution  to  »sz  from 
this  back  course  region,  which  will  be  denoted  by  hs2,  can  be 
written  as  follows: 

J f <*r*>  m .ik<D*R>dydx,  (24) 


where 


d * (x2  ♦ (y-yj2  ♦ h2J1/2. 


R - [ (Xj-x)2  ♦ y2  ♦ zxZ]1/Z.  (26) 

As  before,  the  y- integration  is  performed  by  the  method  of 
stationary  phase  leaving  an  integral  of  the  following  form: 


hsz<V 


-k2hJ0A1/2ei,,/4  f 


* 


Fx(0)  - Cj(0)  lAjfO)  ♦ B^O)] 
f(D(0))x. 

G^O)  « 

D(0)R2(0)S11/2(0) 

-IA1(0)*b1(0)1  x. 

F «(0)  - i i — - (29) 

1 Fj (0)  Aj (0) 

Using  Equation  (28)  to  represent  the  contribution  to  the 
total  scattered  field  from  the  region  behind  the  antenna  and 
applying  Equation  (21)  to  each  illuminated  terrain  segment  lying 
between  the  dipole  and  the  receiver,  we  obtain  the  total  z-cora- 
ponent  of  the  scattered  magnetic  field  by  summing  all  these 
contribution^ : 


H 


sz 


(30) 


The  procedure  for  calculating  the  difference  in  depth  of 
modulation  (DDM)  at  a given  receiver  location  is  quite  straight- 
forward. First,  using  the  procedure  outlined  above,  the  complex 
field  amplitudes  at  the  receiver  due  to  each  element  in  the 
glide  slope  array  are  calculated.  Let  Hc>  Hjjq,  and  denote 
the  total  complex  field  amplitudes  (direct  plus  scattered  summed 
for  all  dipoles  in  the  a-  *ay)  at  the  carrier  frequency,  the  150 
Hz  modulated  frequency,  . 90  Hz  modulated  frequency,  respect- 
ively. In  terms  of  these  complex  amplitudes,  the  difference  in 
depth  of  modulation  is  given  by: 


(31) 


where  Re  denotes  the  real  part  of  the  complex  number.  Full 
scale  deflection  (ISO  microamperes)  corresponds  to  a D.D.M.  of 
.175  for  the  glide  slope.  Consequently,  the  C.D.I.  (course 
deviation  indication)  is  given  by: 


857.14  Re 


C.D.I. 


(32) 


3.  NUMERICAL  RESULTS 

3.1  INTRODUCTION 

The  formalism  developed  in  the  preceding  section  has  been 
applied  to  a whole  series  of  glide  slope  siting  problems.  All 
three  main  types  of  image  glide  slepe  arrays  (null  reference, 
sideband  reference,  capture  effect)  have  been  treated.  Compar- 
ative results  will  be  presented  illustrating  the  performance  of 
each  array  for  a few  terrain  configurations  suggested  to  us  for 
study  by  the  FAA. 

The  null  reference  array  is  the  simplest  of  the  image  glide 
slope  systems.  It  consists  of  two  transmitting  antennas  whose 
heights  are  in  the  ratio  of  2:1.  The  upper  antenna  is  fed  side- 
band only  signal,  the  150  Hz  and  the  90  Hz  components  being  of 
equal  amplitude  but  180°  out  of  phase.  The  lower  antenna  is 
fed  both  carrier  and  sideband.  The  carrier  signal  is  nominally 
40*.  modulated.  The  ratio  of  the  sideband  currents  in  the  upper 
antenna  to  the  sideband  currents  in  the  lower  antenna  is  typically 
0.3.  The  sideband  ratio  of  0.3  would  nominally  yield  3 1.4® 
course  width  (full  deflection  0.7°  above  or  below  the  glide 
path).  Let  Ic,  and  1^^  denote  the  current  amplitudes  fed 

to  the  antennas  at  the  carrier  frequency,  the  150  Hz  modulated 
frequency,  and  the  90  Hz  modulated  frequency,  respectively.  The 
null  reference  array  parameters  discussed  above  can  be  summarized 
as  follows: 


Carrier  Antenna 


Sideband  Antenna 


Height  = h 


I 

I 

I 


c 

150 

90 


1 

» 0.4 
0.4 


Height  = 2h 


I = 0 

r* 

w 

I, 


I 


150 

90 


0.12 

-0.12 


Note  that  all  current  amplitudes  have  been  normalized  relative 
to  the  carrier  current  amplitude  (Ic  = 1). 

The  sideband  reference  array  employs  two  transmitting  dipoles 
whose  heights  are  in  the  ratio  of  3:1.  If  the  lower  antenna  is 
positioned  at  1/2  the  height  of  the  lower  antenna  of  the  null 
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reference  array  and  if  the  upper  antenna  is  positioned  at  3/4 
of  the  height  of  the  upper  antenna  of  the  null  reference  array, 
the  same  glide  angle  is  produced.  Modulated  carrier  (401  modulated) 
is  fed  to  the  lower  antenna.  Both  antennas  are  fed  separate 
sideband  signal.  The  separate  sideband  signals  fed  to  the  two 
antennas  are  equal  in  amplitude  but  are  180°  out  of  phase.  The 
amplitude  ratio  of  the  separate  sideband  signal  to  the  carrier 
sideband  signal  is  typically  0.3.  This  ratio  produces  a nominal 
course  width  of  1.4°  as  in  the  case  of  the  null  reference  array. 

The  sideband  reference  parameters  are  summarized  below. 


Lower  Antenna 


Height  = h 


I 

I 

I 


c 

ISO 

90 


1 

* 0.28 
0 .52 


Upper  Antenna 


Height  * 3h 


11S0 


'90 


0.12 

-0.12 


Note  that  the  I15Q  and  Ig0  current  amplitudes  given  above  for  the 
lower  antenna  represent  the  sums  of  the  carrier  sideband  and 
separate  sideband  signals  U15q  * .4-. 12  * .28,  IgQ  * .4*. 12  = .52). 
All  current  amplitudes  have  been  normalized  relative  to  the  carrier 
signal  amplitude. 

The  capture  effect  glide  slope  array  consists  of  three 
transmi ‘ xing  antennas  whose  heights  are  in  the  ratios  of  1:2:3. 

If  the  lower  andrmiddle  antennas  are  set  at  the  same  heights  as 
the  null  reference  antennas,  the  same  glide  angle  is  produced. 

We  will  not  treat  the  clearance  signal  which  provides  a strong 
fly  up  signal  at  low  approach  angles  but  has  little  effect 
upon  the  glide  angle  and  course  width.  Concerning  the 
primary  signal,  modulated  carrier  is  fed  to  both  the  lower  and 
the  middle  antennas.  The  modulated  carrier  fed  to  the  middle 
antenna  has  half  the  amplitude  and  is  180°  out  of  phase  with  the 
modulated  carrier  fed  to  the  lower  antenna.  The  carrier  signals 
arc  nominally  40*  modulated.  In  addition,  all  three  antennas 
are  fed  separate  sideband  signal.  The  separate  sideband  signals 
fed  to  the  lower  and  upper  antennas  have  half  the  amplitude  and 
are  180°  out  of  phase  with  the  separate  sideband  signal  fed  to 


Note  that  the  I150  and  Igo  current  aaplitudes  given  above  represent 
the  suns  of  the  carrier  sideband  and  separate  sideband  signals. 
Again,  all  values  have  been  normalized  relative  to  the  carrier 
amplitude  in  the  lower  antenna. 


It  should  be  noted  here  that  our  computer  progran  autoaatically 
offsets  the  eleaents  of  each  glide  slope  array  to  correct  for  the 
effects  of  proximity  phase  lag.  The  y coordinates  of  the  eleaents 
are  adjusted  so  that  all  dipoles  are  at  the  saae  slant  distance 
from  the  touchdown  point  on  the  runway  directly  opposite  the  array. 
For  each  array,  one  eleaent  is  held  fixed  while  the  other  eleaents 
are  offset  relative  to  the  fixed  one.  Fcr  the  null  reference  and 
sideband  reference  arrays,  the  lower  element  is  held  fixed  while 
for  the  capture  array,  the  aiddle  eleaent  is  held  fixed.  The 
amount  of  the  offset  for  each  dipole  can  be  calculated  approxim- 
ately in  the  following  manner.  Let  ya  and  h denote,  respectively, 
the  /-coordinate  and  height  of  the  fixed  antenna  element.  The 
/-displacement,  c,  of  any  other  element  in  the  array  relative  to 
the  fixed  one  is  given  approximately  by: 


■^gS  , "v",  " j ’lmt  - 


fc 


c 


h2-H2 


"Zy 


a 


(33) 


where  H is  the  height  of  the  element  which  is  to  be  offset. 

For  the  present  study,  the  wavelength  A was  set  at  3 feet. 

All  arrays  were  positioned  400  feet  from  the  centerline  of  the 
runway.  The  heights  of  the  array  elements  were  set  at  14.33  feet 
and  28.66  feet  (approximately  5A  and  10A)  for  the  null  referer.ee 
array,  7.17  feet  and  21.5  feet  (approximately  2.SA  and  7.5A)  for 
the  sideband  reference  array,  and  14.33  feet,  28.66  feet,  and 
42.99  feet  (approximately  SA,  10A,  and  ISA)  for  the  capture  effect 
array.  If  the  ground  plane  were  perfectly  flat,  a glide  angle  of 
3*  would  be  produced  by  all  three  arrays. 

Two  types  of  data  will  be  presented.  CDI  as  a function  of 
horizontal  distance  from  touchdown  will  be  plotted  for  an  aircraft 
in  level  flight  at  an  altitude  of  1200  feet  as  measured  from  the 
base  of  the  antenna  array  (Zj  * 1200  ft.).  This  horizontal  flight 
data  is  useful  for  studying  course  linearity.  In  addition,  CDI  as 
a function  of  horizontal  distance  from  touchdown  will  be  plotted 
for  an  aircraft  flying  the  nominal  3*  glide  path.  On  the  graphs 
to  be  presented,  positive  CDI  indicates  that  the  150  Hz  modulition 
is  dominant  (fly  up  indication)  while  negative  CDI  indicates  that 
the  90  Hz  modulation  is  dominant  (fly  down  indication). 


3.2  FLAT  TERRAIN 

Figures  2 and  2A  show  the  level  flight  and  nominal  glide 
path  CDI  for  the  null  reference  array  assuming  a perfectly  flat 
ground  plane.  The  data  was  generated  by  summing  the  contributions 
to  the  scattered  field  from  the  first  5000  feet  of  ground  plane. 
The  results  for  all  three  arrays  are  essentially  identical  in  this 
case  so  that  we  will  present  only  the  null  reference  data.  The 
dashed  lines  in  the  figures  represent  the  ideal  results  which  are 
based  upon  simple  image  theory  assuming  a flat,  infinite  ground 
plane.  The  ideal  results  and  the  predictions  of  our  model  are 
clearly  in  close  agreement.  Obviously,  5000  feet  of  flat  ground 
arc  quite  sufficient  to  produce  an  excellent  course.  The  small 
discrepancies  between  the  ideal  results  and  the  predictions  of  our 
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Figure  2.  Null  Reference  Array,  Level  Run  Flight,  Flat  Terrain 


Figure  2A.  Null  Reference  Array, Flyabi 1 i ty  Run, Flat  Terrain 


% 


model  which  are  nost  apparent  in  Figure  2A  result  fron  the  trun- 
cation of  the  ground.  These  deviations  however  amount  to,  at  most, 
approximately  2 nicroaaps.  Figure  2A  shows  a perfect  flare  or 
hyperbola  being  flown  as  would  appear  in  a standard  radio  theodo- 
lite flight  inspection  run.  (Had  a straight  line  glide  path  at  3° 
been  flown,  the  CDI  would  have  assumed  larger  and  larger  positive 
values  (fly  up  indication)  as  the  receiver  approached  the  touch- 
down point  due  to  the  fact  that  the  straight  line  glide  path  is 
actually  only  asymptotic  to  and  actually  below  the  true  zero  DDM 
trajectory  which  is  a hyperbola.) 

3.3  TERRAIN  DISCONTINUITY 

The  next  series  of  graphs  (Figs.  3A,  3B,  SC,  3D,  3E,  3F) 
give  the  predictions  of  our  model  for  a terrain  which  slopes 
downward  away  from  the  array  for  1200  feet  at  an  angle  of  0.6° 
from  the  horizontal  and  then  levels  off  for  an  additional  5S00 
i feet  at  which  point  we  truncate  the  ground  plane  (see  Fig.  3).  The 

dashed  lines  represent  the  ideal  image  theory  predictions  based 
upon  the  infinite,  flat  ground  plane.  The  results  for  all  three 
arrays  (null  reference,  sideband  reference,  and  capture  effect) 
are  very  similar.  The  courses  undergo  a downward  shift  of  0.6° 
relative  to  the  horizontal. 

Referring  to  the  level  run  flight  data  (Figs.  5A,  5B,  3C)  , 
it  will  be  noted  that  the  zero  crossings  occur  at  approximately 
29,500  feet  frGm  touchdown.  For  a receiver  at  a height  of  1200 
feet  as  measured  from  the  base  of  the  antenna  array  (z,  = 1200  ft.), 
the  zero  cross’ng  at  29,500  feet  would  indicate  a glide  angle  of 
approximately  2-35°  relative  to  the  horizontal  (the  x-y  plane). 

The  nominal  glide  angle  is  3°. 

The  flyability  runs  (Figs.  3D,  3E,  3F)  also  reflect  this  down- 
ward shift.  Flying  the  nominal  3°  glide  path,  the  data  shows  a 
nearly  constant  fly  down  indication  of  between  111  and  140  micro- 
amperes. Assuming  a course  width  of  1.4°  (0.7°  above  or  below  the 
glide  path  corresponding  to  full  scale  deflection  of  150  micro- 
amperes) , fly  down  signals  in  this  range  indicate  that  the  receiver 
is  about  0.6®  above  the  true  glide  angle. 
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Figure  3C.  Capture  Kffect  Array,  Level  Run  Flight,  Terrain  Discontinuity 
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Figure  3 1> . Null  Reference*  Array,  Flyability  Run,  Terrain  Discontinuity 
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3.3.1  Adjusted  Antenna  Element  Heights 

In  Oider  to  correct  for  this  downward  shift  of  the  glide  path 
angle,  the  antenna  array  elements  have  been  repositioned  downwards 
to  give  a 3*  glide  path  angle.  This  glide  path  angle  is  determined 
as  the  arithmetical  mean  over  the  distance  between  the  middle 
marker  (3500  feet  from  the  threshold)  and  four  miles  from  the 
threshold  along  the  runway  centerline  extended. 

The  null  reference,  sideband  reference  and  capture  effect 


I 1 

antenna  element 

re-positioning  are  shown  in 

Table  1. 

1 

4 

ft  t 

TABLE  1. 

f ! 

Original 

Positions  (feet) 

New  Positions  (feet) 

j 

■j 

?:  1 
Sf 

Lowest 

Middle  Highest 

Lowest 

Middle  Highest 

* 

[ : 
J.# 

Null  Reference 

14.33 

28.66 

11.99 

23.97 

i 

1 

a 

| 

S.B.  Reference 

7.17 

21.5 

5.88 

17.64 

=; 

% 

Capture  Effect 

14.33 

28.66  42.99 

11.66 

23.32  34.97 

* 

i 

1 

The  level 

flight  results  with  these  new  antenna 

element 

1 

positioning  are  shown  in  Figures  3G,  3H,  and  31.  The  crossovers 
occur  very  nearly  at  the  nominal  glide  angle  of  3°  (which  for  the 
1200  toot  high  horizontal  flight  corresponds  to  approximately 
22,900  feet  on  the  abcissa  of  the  graph). 

The  flyability  run  results  with  this  new  positioning  are  shown 
in  Figures  3J,  3K,  and  iL.  For  all  three  antennas,  the  CDI  has 
been  reduced  to  just  a few  microamperes  by  the  above  repositioning 
of  the  antenna  element  heights. 

3.4  20- FOOT  DRo?0FF 

The  next  series  of  graphs  (Figs.  4A,  4B,  4C,  4D,  4E,  4F)  show 
the  results  for  a terrain  which  extends  flat  for  1200  feet  in 
front  of  the  array  and  then  drops  off  2d  feet  to  a lower  plateau 
which  then  extends  outward  an  additional  3800  feet  (see  Fig.  4). 
From  Figures  4A  and  4D,  it  can  be  seen  that  the  null  reference 
course  tends  to  be  elevated  in.  the  presence  of  such  a terrain. 
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Figure  4A.  Null  Reference  Array,  Level  Run  Flight,  20-Foot  Dropoff 
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Figure  4C . Capture  Effect  Array,  Level  Run  Flight,  20-Foot  Dropoff 
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Figure  4I:.  Capture  Hffect  Array, Fcasihi 1 ity  Run,  20-Foot  Dropoff 


The  horizontal  flight  curve  is  shifted  relative  to  the  ideal  curve 
yielding,  over  much  of  the  course,  stronger  fly  up  signals,  weaker 
fly  down  signals,  and  a zero  crossing  slightly  closer  to  touchdown 
(an  elevated  glide  path)  than  are  predicted  by  image  theory  assum- 
ing a flat  infinite  ground  plane. 

The  nominal  flyability  run  (Fig.  4D)  shows  a consistent  fly 
up  indication  over  most  of  the  course  ranging  between  6 and  13 
microamperes  except  between  the  middle  marker  and  threshold  where 
a fly  down  signal  of  -6  microamperes  occurs. 

The  sideband  reference  data  are  depicted  in  Figures  4B  and  4E. 
Note  in  Figure  4B  that  the  sideband  reference  zero  crossing  point 
is  unaffected  by  the  presence  of  the  dropoff.  Inside  the  zero 
crossing  point,  the  predicted  CDI  is  somewhat  weaker  than  is  pre- 
dicted by  image  theory  while  beyond  the  zero  crossing  point  the 
predicted  CDI  is  somewhat  stronger  than  is  predicted  by  image 
theory,  at  least  in  the  relevant  course  region  (CDI  <150  micro- 
amperes) . 

The  flyability  run  (Fig.  4E)  exhibits  alternate  fly  down  and 
fly  up  indications  with  peak  values  ranging  from  -5  to  *12 
microamperes. 

The  capture  effect  results  are  shown  in  Figures  4C  and  4F. 

The  level  run  flight  exhibits  relatively  close  agreement  between 
the  predicted  CDI  and  the  ideal  image  theory  predictions  except  in 
the  region  beyond  approximately  27,000  feet  where  the  predicted 
fly  up  indication  is  considerably  weaker  than  is  predicted  by 
image  theory. 

The  nominal  flyability  run  exhibits  excursions  ranging  from 
-4  to  7 microamperes.  All  these  deviations  from  ideal  performance 
result  from  reflections  from  the  lower  plateau  which  tend  to 
disrupt  the  nominal  space  modulation  pattern. 

3.4.1  Adjusted  Antenna  Element  Heights 

A 3°  glide  path  angle,  obtained  as  the  arithmetic  mean  over 
the  distance  between  the  middle  marker  and  four  miles  out,  was  next 
flown.  The  3°  glide  path  angle  was  obtained  by  re-positioning 
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the  antenna  element  heights  as  shovm  in  Table  2. 

TABLE  2. 

Original  Positions  (feet)  New  Positions  (feet) 
Lowest  Middle  Highest  Lowest  Middle  Highest 


Null  Reference 

14.33 

28.66 

14.45 

28.89 

S.B.  Reference 

7.17 

21.5 

7.21 

21.62 

Capture  Effect 

14.33  28.66 

42.99 

14.36  28.72 

43.08 

The  level  flight  results  with  this  new  antenna  element  posi- 
tioning are  shown  in  Figures  4G,  4H,  and  41.  These  new  results 
are  only  slightly  different  from  the  original  as  night  be  expected 
with  only  a slight  change  in  the  antenna  element  heights. 

The  flyability  run  results  are  shown  in  Figures  4J,  4K,  and 
4L.  The  null  reference  flyability  run.  Figure  4J  exhibits  a 
mostly  fly  up  indication  with  a peak  value  of  8 microamperes.  Near 
the  middle  marker  region  there  is  a fly  down  indication  with  a peak 
value  of  -11  microamperes. 

The  sideband  reference  and  capture  effect  antennas  (Figures  4K 
and  4L,  respectively)  behave  similarly,  with  peak  values  ranging 
from  -9  to  *8  microamperes  for  the  sideband  antenna  and  ranging 
from  -6  to  *6  microamperes  for  the  capture  effect  antenna. 

3.5  DROPOFF  AND  UPGRADE 

The  next  terrain  configuration  to  be  discussed  is  depicted  in 
Figure  5.  The  terrain  consists  of  1200  feet  of  flat  ground  ending 
in  a 40  foot  drop  followed  by  an  upgrade  rising  58  feet  in  1900 
feet  and  finally  terminating  with  1900  feet  of  flat  ground.  This 
terrain  configuration  as  well  as  the  others  previously  discussed 
was  suggested  to  us  for  study  by  the  FAA.  The  null  reference  data 
are  shown  in  Figures  5A  and  50.  The  level  run  flight  results 
(Fig.  5A)  show  considerable  deviation  from  the  ideal.  The  zero 
crossing  point  is  shifted  by  nearly  2000  feet  toward  touchdown 
(equivalent  to  an  elevation  of  the  glide  path  to  3.19°).  Course 
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Adjusted  Sideband  Reference  Array,  Level  Run  Vllftht,  20-l;oot  Dropoff 
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Figure  4K.  Adjusted  Sideband  Reference  Array,  Flyability  Run, 
20 -Foot  Dropoff 


ed  Capt 


Figure  SB.  Sideband  Reference  Array,  Level  Run  Flight,  Dropoff  and  Upgrade 


Figure  5C.  Capture  Effect  Array,  Level  Run  Flight,  Dropoff  and  Upgrade 


Figure  SI).  Null  Reference  Array,  Flyability  Run,  Dropoff  and  Upgrade 


Sideband  Reference  Array,  Flyability  Run,  Dropoff  and  Upgrade 
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Figure  5F.  Capture  Effect  Array,  Flyahility  Run,  Dropoff  and  Upgrade 


linearity  is  greatly  distorted.  Note  that  the  curve  is  aluost  level 
from  25,000  feet  out  to  29,000  feet. 

The  noainal  flyability  results  (Fig.  5D)  reveal  a strong  fly 
up  indication  of  about  60  aicroaaperes  over  auch  of  the  course  and 
an  even  stronger  fly  down  signal  at  about  7000  feet  froa  touchdown 
of  about  -70  aicroaaperes. 

The  sideband  reference  results  (Figs.  SB  and  5E)  are  qualita- 
tively very  siailar  to  the  null  reference  results  but,  quantita- 
tively, the  distortions  are  not  as  severe.  It  will  be  noted  in 
Figure  5E  that  the  excursions  are  about  *35  aicroaaperes  in  magni- 
tude as  compared  with  60  and  -70  aicroaaperes  in  the  case  of  the 
null  reference  array. 

The  capture  effect  results  (Figs.  SC  and  SF)  show  the  least 
distortion.  In  Figure  SC,  it  can  be  seen  that  the  z ;ro  crossing 
point  undergoes  a relatively  aodest  displacement , 'he  noainal  glide 
path  results  (Fig.  SE)  exhibits  excursions  ranging  from  -11  to  17 
microamperes,  a decided  iaproveaent  over  the  null  reference  and 
sideband  reference  arrays. 

3.S.1  Adjusted  Antenna  Element  Heights 

In  an  attempt  to  improve  the  flyability  of  this  course,  the 
antenna  element  heights  were  adjusted  to  give  a mean  3°  glide  path 
angle  over  the  distance  between  the  middle  marker  and  four  miles 
out.  The  new  antenna  element  heights  arc  shown  in  Table  5 

TABLE  3. 


Original 

Positions  (feet) 

New  Positions 

( feet) 

Lowest 

Middle 

Highest 

Lowest 

Middle 

Highest 

Null  Reference 

14. -<3 

28.66 

IS. 02 

50.04 

S.B.  Reference 

7.17 

21. S 

7.35 

22.04 

Capture  Effect 

14.33 

28.66 

42.99 

14. S9 

29.18 

45.77 

The  level  flight  results  with  this  new  positioning  arc  shown  in 
Figures  5G,  5H,  and  SI.  The  flyability  results  are  shown  in 
Figures  5J,  5K,  and  5L. 
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Figure  50.  Adjusted  Null  Reference  Arruy,  Level  Run  Flight,  Dropoff  and  Upgrade 
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PJ HUM*  5M.  Adjust fd  Sideband  Refereme  Array,  Level  Run  l;  1 i >\l»t  . Dropoff  and  Upgrade 
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Figure  Sk . Adjusted  Sideband  Reference  Array,  FlyabUlty  Run,  Dropoff  and  Upgrade 
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Re-positioning  of  the  antenna  element  heights  for  this  type  of 
terrain  to  obtain  a mean  3°  glide  path  angle  does  not  improve  the 
course  though  the  cross  overs  (Figures  5G,  5H,  and  51)  do  occur 
somewhat  closer  to  the  required  3°  glide  angle.  As  can  be  seen  on 
comparing  Figure  5J,  5K,  and  5L  with  the  uncorrected  values  shown 
in  Figures  Sb  5E,  and  5F,  the  flyability  of  the  course  has  not 
improved  by  this  re-positioning.  The  peak-to-peak  values  are  now 
-115  to  21  microamperes  for  the  null  reference  antenna;  -57  to  21 
microamperes  for  the  sideband  reference  antenna;  and  -22  to  7 micro- 
amperes for  the  capture  effect  antenna.  While  the  fly  up  indica- 
tions have  in  all  cases  been  reduced,  the-  fly  down  indications  have 
all  gone  up.  This  is  because  the  re-positioning  of  the  antenna 
element  heights  simply  shifts  the  CDI  pattern  either  up  or  down. 

This  type  of  correction  is  very  useful  when  a uniform  shift  of  the 
CDI  pattern  occurs,  for  example,  when  the  terrain  has  an  up  or  down 
slope,  as  v;,s  the  case  depicted  in  Figure  3.  Adjusting  the  antenna 
element  heights  there  for  a mean  3°  glide  path  angle  shifted  the 
whole  pattern  back  toward  zero  CDI.  However,  for  the  terrain 
depicted  in  Figure  5,  which  produced  both  up  and  downward  CDI 
shifts  (as  shown  in  Figures  5D,  5E,  and  5F) , tne  simple  translation 
of  the  CDI  pattern  that  was  produced  by  adjusting  the  antenna 
element  heights  for  a mean  3°  glide  path  angle,  cannot  yield  zero 
CDI  over  the  whole  course.  By  reducing  the  fly  up  indications, 
the  fly  down  indications,  in  fact,  became  stronger. 


We  do  emphasize  again,  however,  that  for  this  terrain  the 
sideband  reference  antenna  produced  a much  improved  course  over 
the  null  reference,  and  that  the  capture  effect  antenna  a much 
improved  course  over  the  siaerand  reference  ante.ma. 


3.6  UPGRADE,  DROPOFF  AND  UPGRADE 

The  next  terrain  configuration  to  be  discussed  is  depicted  in 
Figure  6.  The  terrain  consists  of  1000  feet  of  ground  sloping  up 
at  3/4°  ending  in  a 20  foot  drop  which  is  then  followed  by  500  feet 
of  flat  ground  followed  by  3500  feet  of  ground  again  sloping  up  at 
3/4°.  The  results  are  shown  in  Figures  6A,  6B,  6C,  6D,  6E,  and  6F. 
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6.  Upgrade, 


Null  Reference  Array.  Flyability  Run,  Upgrade,  Dropoff  and  Upgra 
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Figure  6F..  Sideband  Reference  Array,  Flyability  Run,  Upgrade,  Dropoff  and  Upgrade 


Figure  6F . Capture  liffcct  Array,  Flyability  Run,  Upgrade,  Dropoff  and  Upgrade 


Looking  at  the  level  flight  runs  first.  Figures  6A,  6B»  and 
6C,  we  see  that  the  null  reference  antenna  data  exhibits  no  zero 
cross  over  in  the  range  considered  (though  there  is  a probable  zero 
cross  over  at  an  elevated  glide  path  angle  of  3.7°;  in  contrast  to 
this  the  sideband  reference  antenna  has  two  zero  cross  overs  in 
the  same  range,  one  at  about  2.1°  and  the  other  at  approximately 
2.5°  (and  another  probable  zero  crossover  at  about  3.7°)  to  pro- 
duce a highly  non  linear  distorted  course;  the  capture  effect 
antenna  appears  to  have  two  zero  cross  overs  one  at  1.5°  and  an- 
other at  3.6®. 

The  flyability  results  are  shown  in  Figures  6D,  6E,  and  6F. 

The  null  reference  data  (Figure  6D)  reveals  a strong  fly  up  signal 
of  about  15S  microamperes  which  is  equivalent  to  almost  3/4®,  which 
is  just  the  angle  at  which  the  two  separate  pieces  of  ground  slope 
upwards  (see  Figure  6).  Clearly,  to  the  null  reference  antenna, 
this  typr*  of  terrain  appears  as  a simple  upslope. 

Interestingly,  the  sideband  reference  results  (Figure  6E)  for 
this  type  of  terrain  are  even  pooler  than  the  preceding  for  the 
null  reference  antenna. 

The  best  performance  is  exhibited  by  the  capture  effect 
antenna  whose  results  are  depicted  in  Figure  6F.  There  is  a fly 
up  signal  over  the  flight  path  whose  peak-to-peak  excursions  range 
from  46  microamperes  to  90  microamperes,  about  half  that  of  the 
null  reference  antenna.  The  terrain  does  not  appear  as  a simple 
3/4°  upslope  to  the  capture  effect  antenna  which  would  have  then 
produced  the  same  155  aicroaap*-  • • excursion  as  it  did  for  the  null 
reference  antenna  (see  for  exaa.>ie,  the  results  for  the  simple 
upslope  pictured  in  Figure  3). 

3.6.1  Adjusted  Antenna  Element  Heights 

Figures  6G,  6H,  61,  6J,  6K,  and  6L  exhibit  the  results  for 
the  terrain  of  Figure  6 after  the  antenna  element  heights  have  been 
adjusted  to  provide  a 3®  glide  path.  This  3®  glide  path  angle  was 
determined  as  the  arithmetic  mean  over  the  distance  between  the 
middle  marker  and  four  miles  from  the  threshold.  The  adjusted 
heights  of  the  elements  are  given  in  Table  4. 
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F:igure  6G.  Ail  justed  Null  Reference  Array,  Level  Run  Plight,  Upgrade,  Dropoff  and  Upgrade 


j*urc  61.  Adjusted  Capture  Kffcct  Array,  l.evel  Run  Plight,  Upgrade, 
Dropoff  and  Upgrade 
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TABLE  4. 

Original  Positions  (feet" 
Lowest  Middle  Highest 


New  Positions  (feet) 
Lowest  Middle  Highest 


The  level  flight  results  with  these  re-adjusted  antenna 
element  heights  are  shown  in  Figures  6G,  6H,  and  61.  These  level 
flight  results  still  show  considerable  deviation  from  the  ideal; 
with  course  linearity,  in  all  three  cases,  distorted,  though  the 
null  reference  antenna  data  (Figure  6G) , does  exhibit  the  correct 
zero  cross  over  at  3° . 

The  sideband  reference  level  flight  results  (Figure  6H) 
exhibits  no  less  than  three  zero  cross  over  points  within  the 
range  considered  here  (20.000*  to  46,000*).  Two  of  these  zero 
crossings  (one  at  about  3°  and  the  other  at  about  2°)  are  in  the 
right  direction,  which  therefore  could  present  a serious  problem 
to  the  pilot.  These  resvilts  indicate  tr.-al  it  is  possible  to  fly 
a zero  CDI  course  at  the  incorrect  glide  angle  using  a sideband 
reference  antenna  over  this  terrain. 

The  level  flignt  values  produced  by  the  capture  effect  system 
are  shown  in  Figure  6' . The  zero  crossing  here  is  at  approximately 
2.8". 

The  flyability  of  rhe  course  has  been  dramatically  improved 
in  all  three  cases  with  the  new  re-positioned  antenna  element 
heights.  The  null  reference  antenna  results,  shown  in  Figure  6J , 
has  excursions  between  * 2C  microamperes  and  +9  microamperes  (com- 
pared to  155  microamperes  before  re-positioning).  This  is  con- 
siderably better  than  the  sideband  reference  results  (Figure  6K) 
which  could  not  be  corrected  to  give  less  than  peak-to-peak  micro- 
ampere excursions  of  -73  to  +25,  and  even  somewhat  better  than  the 
capture  effect  antenna  results  (Figure  6L)  which  show  excursions 


Null  Reference 

14.33 

28.66 

19.28 

38.55  ] 

1 

S.B.  Reference 

7.17 

21.5 

9.87 

29.61  ; 

J 

Capture  Effect 

14.33  28.66 

42.99 

16.48  32.97 

49.45  ! 

of  between  -19  and  +25  microamperes  after  re-positioning  of  the 
antenna  element  heights. 

It  is  interesting  that  for  this  type  of  terrain,  consisting 
of  two  upslopes  and  a discontinuity,  it  has  still  been  possible 
to  correct  the  null  reference  antenna  to  produce  results  at  least 
as  good  as  those  produced  by  the  capture  effect  antenna.  This 
will  be  the  case,  in  general,  whenever  the  terrain  consists  of 
simple  up  or  downslopes. 

3.7  HILL  TERRAIN 

The  final  terrain  configuration  to  be  discussed  is  depicted 
in  Figure  7.  The  terrain  consists  of  3000  feet  of  flat  terrain 
followed  by  an  upgrade  of  5°  extending  for  another  2000  feet.  The 
null  reference  data  are  shown  in  Figures  7A  and  7D.  The  level 
flight  results  (Figure  7A)  show  a zero  crossing  point  which  is 
shifted  about  1000  feet  toward  touchdown,  equivalent  to  an  eleva- 
tion of  che  glidepath  to  3.12°. 

The  flyability  results  (Figure  7D)  reveal  a strong  fly  up 
indication  of  67  microamperes  over  the  course  from  30,000  feet  tr 
about  9000  feet  from  touchdown.  This  is  followed  by  several  very 
strong  fly  down  indications,  the  highest  of  which  occurs  approxi- 
mately 5S00  feet  from  touchdown  with  a peak  value  of  -182  micro- 
amperes. 

The  sideband  reference  results  are  shown  in  Figures  78  and  7E. 
The  level  flight  data  (Figure  7B)  shows  a cross  over  precisely  at 
the  3°  glide  path  angle,  though  an  additional  false  cross  over 
occurs  at  1.6°  (about  43,000  feet  from  touchdown  for  the  1200  foot 
high  level  flight).  The  flyability  data  (Figure  7E)  reveals  an 
almost  distortion- free  course  until  about  18,000  feet  from  touch- 
down. However,  as  the  aircraft  approaches  closer,  a gradually- 
increasing  flv  up  signal  is  received  which  peaks  at  87  microamperes 
at  about  8SOO  feet  from  touchdown.  This  is  followed  by  even 
stronger  fly  down  signals,  the  largest  of  which  peaks  at  -16S 
microamperes  about  5500  feet  from  touchdown. 
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The  capture  effect  antenna  system  results  are  shown  in  Figures 
| 7C  and  7F.  The  level  flight  data  (Figure  7C)  exhibits  a zero  cross 
over  which  is  very  nearly  at  the  correct  3°  glide  path  angle, 
though  here  too,  as  in  the  sideband  reference  data,  there  appears 
to  be  a false  cross  over  (at  about  1.46°).  The  nominal  glide  path 
■ results  (Figure  7F)  show  that  the  capture  effect  antenna  system 

offers  a decided  improvement  over  the  two  previous  antenna  systems, 
i There  is  a near  zero  CDI  over  most  of  the  course,  never  exceeding 

) 27  microamperes,  from  30,000  feet  until  8000  feet  from  touchdown. 

I This  mild  fly  up  signal  is  followed  by  several  fly  down  signals 

! which ^however,  never  exceed  -45  microamperes  (which  is  only  about 

i a quarter  that  of  the  other  two  systems).  Clearly,  for  this  type 

of  terrain,  characterized  by  flat  ground  followed  by  a mildly 

i 

I steep  hill,  only  the  capture  effect  system  can  produce  an  acceptable 

i course  structure.  This  is  true  even  aftei  the  antenna  element 

heights  are  adjusted  to  form  a mean  3°  glide  path  angle. 

3.7.1  Adjusted  Antenna  Element  Heights 

This  mean  glide  path  angle  was  formed  by  re-positioning  the 
antenna  element  heights  as  shown  in  Table  5: 


TABLE  5. 

Original  Positions  (feet)  New  Positions  (feet) 
Lowest  Middle  Highest  Lowest  Middle  Highest 


Null  Reference 

14.33 

28.66 

14.64 

29.28 

S.8,  Reference 

7.17 

21.5 

7.22 

21.7 

Capture  Effect 

14.33  28.66 

42.99 

14.35  28.71 

4 3.07 

The  results  with  this  re-positioning  of  the  antenna  element 
heights  are  shown  ifi  Figures  7G,  7H,  71,  7J,  7K,  and  7L.  The  level 
flight  data,  Figures  7G,  ?H,  and  71,  arc  similar  to  the  previous 
level  flight  data,  which  was  to  be  expected  with  the  modest  change 
in  the  rntenna  element  heights  needed  to  obtain  a mean  3°  glide 
path  angle. 


j 


1 
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Figure  7K . Adjusted  Sideband  Reference  Array,  Flyability  Run,  Hill 


Figure  7L.  Adjusted  Capture  Bffect  Array,  Flyability  Run,  Hill 


The  f lyability  data  are  shown  in  Figures  7J,  7K,  and  7L.  With 
the  null  reference  antenna  (Figure  7J)  there  is  less  of  a fly  up 
indication  after  antenna  adjustment  then  before  (an  almost  unifora 
S6  aicroaaperes  as  coapared  to  6/  aicroaaperes  before  the  re- 
positioning) but  larger  fly  down  peaks  (-208  aicroaaperes  as  con- 
pared  to  -182  aicroaaperes).  Similarly,  for  the  sideband  reference 
antenna:  the  fly  up  signal  decreased  froa  a maximum  of  87  micro- 

amperes to  a aaxiaua  of  82  aicroaaperes,  but  the  fly  down  signal 
increased  from  -165  aicroaaperes  to  -171  aicroaaperes  aaxiaua. 

The  course  structure  using  the  capture  effect  system  remains 
essentially  the  saae  as  before  re-uositioning. 

As  also  found  for  the  case  of  the  terrain  depict in  Figure 
5,  a simple  re-positioning  of  the  antenna  element  heights  to 
obtain  a aean  3°  glide  path  angle  over  the  distance  between  the 
aiddle  marker  and  four  ailes  out  does  not  iaprove  the  course 
structure  for  the  terrain  of  Figure  7.  As  before,  here  too,  the 
original  course  structure  showed  both  fly  up  and  fly  down  excur- 
sions. In  correcting  the  one,  the  other  becomes  worse.  The  only 
system  capable  of  producing  an  acceptable  course  structure  for 
this  ’‘hill"  type  terrain  is  the  capture  effect  system. 
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PART  II.  USER'S  MANUAL* 

1.  introduction:  definition  of  instrument  landing  system 


The  ILSGLD  program  has  been  written  to  simulate  certain 
terrain  conditions  which  affect  the  glide  slope  portion  of  the 
Instrument  Landing  System.  The  ILS  is  used  to  provide  signals  for 
the  safe  navigation  of  landing  aircraft  during  periods  of  low 
cloud  cover  and  other  conditions  of  restricted  visual  range. 
Separate  systems  are  used  to  communicate  vertical  and  horizontal 
information;  the  vertical  system  is  called  the  "glide  slope". 

This  system  operates  by  the  transmission  of  an  RF  carrier, 
amplitude  modulated  by  two  audio  frequencies,  beamed  to  approach- 
ing airborne  receivers.  In  an  instrumented  aircraft,  the  glide 
slope  receiver  serves  to  demodulate  the  RF  signal,  amplify  and 
isolate  the  corresponding  audio  signals  and  derive  a signal  to 
drive  the  ILS  vertical  display  in  the  cockpit.  The  pilot,  by 
reading  the  display,  can  determine  if  he  is  on  course,  above  or 
below  the  glide  path.  These  signals  must  be  strong  enough  to 
cover  a radius  of  15  miles  around  the  antenna. 

The  directional  information  is  determined  by  the  relative 
strengths  of  the  transmitted  sideband  signals.  The  audio  fre- 
quency modulations,  which  are  fixed  at  90  Hz  and  150  Hz,  are 
radiated  in  different  angular  patterns  with  respect  to  the  intended 
glidepath.  The  "course"  is  defined  as  the  locus  of  points  where 
the  amplitudes  of  the  two  modulations  are  equal.  The  display  of  a 
difference  of  the  amplitudes  (90  Hz  and  150  Hz)  of  the  sidebands 
is  referred  to  as  the  Course  Deviation  Indication.  Thus,  the  CDI 
is  the  pilot’s  indication  as  to  what  his  altitude  is  relative  to 
the  glidepath.  The  CDI  is  measured  in  microamps.  The  actual 
course  generated  by  any  particular  ILS  installation  will  deviate 
from  the  ideal  due  to  irregularities  in  the  terrain.  The  devia- 
tion of  the  CDI,  caused  by  these  irregularities,  from  t.ie  ideal 
receiver  reading  at  that  point  in  space  (e.g.,  on  the  glidepath 
a CDI  reading  other  than  0)  is  the  derogation  effect. 

•Part  H is  in  actuality  Section  4 of  the  report;  accordingly. 
Chapters  1-5,  Part  II,  are  paginated  as  4-1  through  4-16. 
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The  glide  slope  system  transmits  an  asymmetrical  pattern  by 
beaming  a "carrier  plus  sideband"  pattern  and  a "sideband  only" 
pattern,  the  composite  of  which  gives  the  desired  effect. 
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2.  ANTENNA  PATTERNS 


The  proper  angular  variation  of  the  transaitted  90  Hz  and  the 
150  Hz  modulation  is  achieved  by  the  radiation  of  the  two  independ- 
ent sideband  patterns  by  the  transmitting  antenne*.  One  pattern, 
the  "carrier  plus  sidebands"  (OS)  signal,  is  radiated  in  a 
symmetrical  pattern;  the  other  pattern,  the  "sidebands-only"  (SO) 
signal,  is  radiated  in  an  "anti -sy metrical"  pattern  relative  to 
the  prescribed  glide  angle  (see  Figure  8.) 

The  OS  signal  is  composed  of  a carrier  wave  and  paired  wide- 
band waves  at  90  Hz  and  150  Hz.  The  sideband  amplitudes  are  equal 
and  represent  a 40$  modulation  of  the  carrier  wave  (or  a "depth 
of  modulation"  of  0.4)  at  both  frequencies.  The  SO  wave  is  a 
carrier  wave  that  is  equally  modulated  at  90  Hz  and  ISO  Hz  to  the 
extent  that  it  retains  no  pure  carrier  component. 

The  spatial  modulation  pattern  obtained  by  combining  the 
symmetrical  OS  pattern  with  the  "anti -symmetrical"  SO  pattern  is 
illustrated  in  Figure  8.  At  a given  receiver  point  the  total 
signal  is  the  OS  carrier  plus  the  combined  sideband  amplitudes  of 
the  OS  and  SO  patterns.  The  sideband  amplitudes  are  phased  so 
that  above  the  glide  path  the  90  Hz  amplitudes  add  and  the  150  Hz 
amplitudes  subtract,  while  below  the  glide  path,  the  90  Hz  ampli- 
tudes subtract  and  the  150  Hz  amplitudes  add. 

Any  angular  deviation  of  the  airplane's  receiver  from  the 
correct  course  results  in  a "Difference  in  Depth  of  Modulation" 
(DDM)  between  the  150  Hz  and  90  Hz  signals.  Since  the  strength  of 
OS  and  SO  signals  fall  off  at  the  same  rate  with  distance  from 
the  transmitting  antenna,  the  DDM  is  independent  of  range. 

The  antenna  patterns  are  generated  by  using  arrays  of  one  or 
more  elements  in  combination  with  the  ground  as  reflector.  The 
effect  of  an  ideal  ground  plane  on  a single  element  is  3$  though 
there  was  an  "image"  element  located  below  the  ground  and  radiat- 
ing equal  power  to  the  real  clement  but  with  opposite  phase.  In 
a glide  slope  array  there  will  be  two  or  more  elements  radiating 
different  signals  to  give  the  desired  combined  antenna  pattern. 
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(Note  that  this  means  that  a glide  slope  array  is  defined  by  giving, 
for  each  element,  its  location  and  complex  amplitude'  of  the 
carrier  and  sidebands.)  However,  the  real  ground  is  not  an  ideal 
plane.  This  has  the  effect  of  distorting  the  element  patterns  and 
results  in  a derrogation  of  the  glide  slope  system  performance. 


3.  SIMULATION  DESCRIPTION  OF  1LSGLD 


The  ILSGLD  simulation  program  allows  the  user  to  study  the 
effects  of  various  ground  shapes  and  antenna  types.  Using  this 
program,  the  user  may  determine  what  antenna  systems  will  give 
acceptable  performance  with  a particular  giound  configuration  or 
what  ground  shapes  cause  problems  with  various  antennae. 

The  program  uses  a ground  description,  an  antenna  description, 
jnd  a set  of  spatial  coordinates  of  the  receiver  location  as  input. 
The  program  calculates  and  outputs  for  each  receiver  location,  the 
CDI  at  that  point  both  from  the  "real"  ground  and  from  an  " ideal" 
(perfectly  conducting  and  flat)  ground  plane.  In  addition,  the 
effects  of  the  electrical  and  mechanical  inertia  are  included  and 
"dynamic  CDI's"  are  also  output  for  both  esses. 

The  inputs  for  ILSGLD  are  data  files  prepared  by  the  program 
FMAKE.  The  outputs  from  ILSGLD  are  in  data  files  which  are  used 
as  input  to  a third  program  GLDPLT.  GLDPLT  produces  graphical 
outputs  of  t%.e  various  CDI's.  The  description  and  usage  of  these 
programs  are  described  in  succeeding  sections. 

The  ILSGLD  simulation  makes  certain  simplifying  assumptions. 
These  include: 

a.  perfectly  reflecting  ground  surface 

b.  far  field  scattering  - all  scattering  from  points  on  the 
ground  surface  is  assumed  independent  of  all  other  points, 
thus  multiple  reflections  and  near  field  interactions  arc 
ignored 

c.  noise  free  environment 

d.  relative  field  strengths  - the  absolute  field  strengths 
involved  are  not  calculated.  Thus,  while  the  CDI’s  can 
be  calculated  in  microamperes,  the  absolute  electric 
field  intensities  are  not  ascertained 

e.  idealized  receiver  model 

f.  geometrical  shadowing  * the  shadowing  of  one  portion  of 
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the  ground  on  another  is  done  by  straight  ray  approxima- 
tions assuming  a total  cutoff  with  no  diffraction. 

g.  antenna  elements  are  assumed  to  be  simple  dipoles. 

3.1  METHOD  OF  SIMULATION 

An  antenna  clement  is  described  by  giving  its  x- , y-  and  z- 
conrdinates  and  the  complex  amplitudes  of  its  radiated  amplitude 
at  the  three  frequencies  (carrier,  150  Hz,  and  90  Hz  sidebands). 

The  ground  is  broken  up  into  strips;  these  strips  have  an  infinite 
width  and  a finite  length.  The  infinite  extent  is  parallel  to  the 
y-axis,  that  is  at  a right  angle  to  the  runway  centerline.  Thus,  a 
ground  strip  is  described  by  giving  the  x-  and  z-coordinates  of 
the  leading  and  trailing  edges  of  the  strip.  A receiver  location 
is  described  by  giving  its  time,  x* , y-  and  z-coordinates. 

The  basic  part  of  the  simulation  consists  of  calculating  the 
field  at  a receiver  location.  This  field  is  caused  by  the  power 
radiated  from  an  element  and  reflected  fro*  a strip.  A detailed 
description  of  the  method  of  calculation  has  been  given  previously 
(Part  I).  The  receiver  field  can  be  expressed  as  a complex  gain 
factor  tines  the  radiated  antenna  field.  This  gain  factor  is 
expressed  as  a double  integral  over  the  strip.  The  integration 
along  the  infinite  extent  (y-axis  direction)  was  approximated  by 
using  the  stationary  phase  method.  The  resulting  single  Integra? 
is  solved  by  a modified  trapezoid  rule.  The  trapezoid  rule  is 
used  with  the  spacing  between  sample  points  adjusted  for  the 
derivative  of  the  integrand. 

Thus,  for  a given  receiver  location,  the  program  takes  an 
antenna  element  and  calculates  the  "gain  factor"  for  each  ground 
strip,  multiplies  by  the  radiated  power,  and  sum*  the  three  complex 
field  intensities  for  ail  ground  strips.  The  direct  field  of  the 
element  at  the  receiver  is  then  added.  This  process  is  repeated 
for  all  elements,  giving  the  total  field  at  the  receiver.  Using 
an  idealized  receiver  model  the  voltages  at  the  output  of  the 
receiver  are  calculated  and  the  CD!  derived.  This  Is  then  repeated 
for  the  next  receiver  location  by  summing  the  fields  over  all  the 


elements  and  over  all  the  strips.  For  comparison  purposes,  the 
"ideal"  ground  plane  field  is  also  calculated  and  summed  separately. 
The  CDI's  due  to  actual  and  idealized  ground  conditions  may  then 
be  compared.  Finally,  a second  pair  of  CDI’s  are  derived  by  using 
a first  order  difference  equation  with  the  instrument  time  con- 
stant and  the  time  differc  .ce  from  the  last  point,  thus  simulating 
the  "inertia"  of  the  ILS  glide  slope  receiver  and  display  (so 
called  "dynamic"  CDI). 


At  each  receiver  point,  an  output  record  is  generated  in 
file  STRIP.DAT  consisting  of  the  x- , y-,  z-,  and  time  coordinates 
of  the  receive  plus  the  four  CDI's. 

If  the  terrain  is  sufficiently  irregular,  part  of  the  energy 
radiated  toward  a point  on  the  ground  may  be  intercepted  by 
another  piece  of  the  ground  closer  to  the  antenna.  This  shadowing 
is  complex  including,  as  it  does,  diffraction  as  well  as  reflec- 
tion. This  is  approximated  in  the  simulation  by  using  ray  optics; 
that  is  . the  shadow  is  assumed  to  have  no  diffraction  at  the  edges 
and  a lero  field  amplitude  inside  the  shadow.  The  program  does 
this  by  assuming  the  ground  is  continuous,  i.e.  the  far  edge  of 
one  strip  is  the  near  edge  of  the  next,  and  keeping  track  of  the 
"furthest"  (in  angular  sense)  edge.  If  part  (or  all)  of  the  next 
strip  is  "below"  that  edge  then  that  part  (or  all)  will  not  be 
included  in  the  trapezoid  integration.  For  example  in  the  sketch 
below 


all  of  strip  I,  none  of  strip  II  and  that  part  of  strip  III  between 
D and  E will  be  included  in  the  integration. 

For  the  receiver  antenna  a "semi"  directional  antenna  is 
assumed;  that  is,  only  the  incident  fields  from  the  front  half- 
sphere around  the  receiver  are  included.  This  is  as  though  an 
omni-directional  antenna  was  used  but  blocked  by  the  fuselage  from 
receiving  signals  from  the  direction  of  the  tail.  This  is  done  in 
the  program  by  stopping  the  summation  of  fields  over  the  ground 
strips  at  a point  directly  below  the  receiver. 

The  back  half-plane  is  assumed  to  be  an  ideal  flat  horizontal 
reflector  of  infinite  extent.  The  field  from  this  is  included  by 
adding  the  "gain  factor"  for  this  as  an  initial  strip  to  that  cal- 
culated by  integrating  over  the  "real"  strips. 

5.2  OPERATION 

ILSGLD  assumes  the  ground  description  is  a file  called 
GRND.DAT,  and  that  the  receiver  locations  are  in  a file  called 
PATH.DAT.  The  user  starts  the  program  and  then  inputs  the  name  of 
the  file  containing  the  antenna  description.  The  simulation  will 
be  run,  and  the  output  will  be  found  in  file  STRIP.DAT. 

5.5  CONCLUSION 

ILSGLD  represents  an  initial  effort  at  glide  slope  landing 
system  simulation.  It  has  been  exercised  with  simple  test  cases. 
Two  possible  areas  for  future  work  include: 

1.  terrain  variation  perpendicular  to  runway  centerline 
direction  - the  ground  does  not  usually  have  a constant 
profile  along  the  y-direction  as  is  presently  assumed; 

2.  varying  conductivity  - the  physical  ground  surface  is 
not  necessarily  the  "surface"  as  seen  by  the  radio 
frequency  fields  as  is  presently  assumed. 
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4.  FMAKE  PROGRAM  DESCRIPTION 


FMAKE  is  a file  generation  program  used  to  create  input  files 
for  ILSGLD.  It  is  designed  to  be  used  interactively.  The  user 
starts  by  running  the  program.  The  program  will  respond  by  typing: 

INPUT  SWITCH: 

The  user  then  types  in  a single  character  switch,  followed  by  a 
<CR>,  for  the  file  he  wishes  to  generate.  The  program  will  then 
respond  with  a request  for  the  input  required  for  that  file.  If 
a blank  is  used  as  the  input  switch  the  program  will  terminate. 

If  any  character  other  than  the  ones  explained  below  are  used  an 
error  message  will  be  given.  After  each  file  is  generated,  the 
program  will  return  to  the  switch  input  point  thus  allowing  the 
user  to  generate  the  data  files  for  many  simulation  runs  in  one 
sitting. 

(N.B.  All  units  are  in  feet  unless  otherwise  stated) 

4.1  SWITCH:  Y 

The  program  will  type: 

INPUT  YO,  LAMBDA: 

The  user  then  inputs  in  free  field  format  the  y-offset  (i.e.  the 
y-coordinate)  of  the  antenna  elements  and  the  wave  length  of  the 
carrier,  followed  by  a <CR>.  The  y-offset  is  the  distance  from 
the  base  of  the  antenna  to  the  centerline  of  the  runway.  As  this 
information  is  required  for  both  the  antenna  description  and  the 
flight  path  it  is  input  with  a separate  switch  to  avoid  repetition. 

4.2  SWITCH:  G 

This  is  used  to  input  the  ground  description.  The  program 
will  respond  with: 

INPUT  CROUND  FILE  NAME: 

The  user  then  inputs  a <five(S)  character  name  for  the  ground 
descriotion  file,  followed  by  a <CR>.  ILSGLD  requires  the  ground 
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file  to  be  called  GRND.DAT  (The  .DAT  extension  is  a system 
default).  FMAKE  allows  the  user  to  generate  several  ground  files 
with  different  names  at  one  sitting.  Then  by  using  system  renaming 
commands  a single  batch  job  may  be  set  up  that  will  run  many 
simulations  without  further  user  interactions. 

The  program  will  then  type: 

INPUT  GROUND  LABEL 

The  user  then  inputs  up  to  40  characters  to  be  used  as  a label  for 
this  ground  description.  This  label  is  carried  in  the  ground 
description  file  and  will  be  placed  in  the  output  file  by  ILSGLD. 
This  allows  the  GLDPLT  program  to  label  the  plots  with  the  ground 
description.  This  is  necessary  if  a batch  job  generates  plots  from 
more  than  one  simulation. 


The  program  will  then  type: 

INPUT  GROUND  SEGMENTS.  STARTING  FROM  ANTENNA,  GIVE 
CONSECUTIVELY  EITHER  X AND  Z INCREMENTS,  OR  THE  LENGTH  AND 
ANGLE  FROM  HORIZONTAL  IN  DEGREES,  SEPARATED  BY  A ZERO.  HIT 
CARRIAGE  RETURN  FOR  END  OF  DATA,  OR  IF  THERE  ARE  NO  MORE 
SrRIPS. 

j'he  user  inputs  ground  strip  end  points  in  either  Cartesian  or 
polar  increments  using  two  or  three  fields  respectively.  These 
are  input  in  floating  free  field  format  followed  by  a <CR>.  The 
first  point  is  taken  relative  to  the  origin  and  is  the  near  edge 
of  the  closest  strip  to  the  antenna.  The  second  point  is  the  far 
edge  of  the  first  strip  and  the  near  edge  of  the  second  strip. 
This  second  point  is  relative  to  the  first.  This  will  continue 
for  additional  strips  until  a <CR>  is  hit  with  no  data,  at  which 
point  the  program  will  return  to  the  switch  input.  For  example 
to  input  this  profile: 
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The  input  would  be 


0.  ,0. 

1200.  , 0. , - .6 

0.,  *100. 

2400.,  170. 

There  is  a maximum  of  20  strips  allowed  in  both  FMAKE  and 
ILSGLD.  This  is  determined  by  array  sizes  and  could  be  changed  if 
desired. 

4.3  SWITCH:  P 

This  is  used  to  generate  a flight  path  file.  The  program 
will  respond: 

INPUT  FLIGHT  PATH  FILE  NAME: 

ihe  user  then  inputs  a five(5)  character  file  name  (the  name  must 
be  exactly  S characters).  ILSGLD  requires  the  flight  path  to  be 
in  a file  called  PATH . DAT  (for  explanation  of  multiple  files  see 
SWITCH :G).  The  program  will  then  type: 

INPUT  FLIGHT  PATH  TITLE: 

The  user  then  inputs  a title,  using  up  to  40  characters,  for  the 
flight  path.  This  is  used  as  a label  on  the  plots  output  by 
GLDPLT.  The  program  will  then  type: 

INPUT  FLIGHT  PATH  TYPE: 

The  user  has  a choice  of  two  flight  path  types,  linear  or  hyperbolic. 
For  a linear  flight  path,  type  a <CR>;  for  a hyperbolic  path  type 
a *G'  followed  by  a <CR>.  When  it  is  a linear  flight  the  program 
will  respond: 

INPUT  X0,  Y0,  20: 

The  user  then  inputs  the  x- , y- , and  z-coordinates  of  the  first 
receiver  point  in  free  field  floating  point  format  followed  by  a 
<CR>.  The  program  will  respond: 

INPUT  XF,  YF,  2F : 
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The  user  then  inputs  the  final  receiver  location  in  the  same  way. 
The  program  then  types: 

INPUT  * OF  POINTS,  VELOCITY.  TIME  CONSTANT: 


The  user  then  inputs  the  total  number  of  receiver  locations 
desired  in  integer  free  field  format,  followed  by  the  velocity 
of  the  aircraft  in  feet/sec.  in  floating  point  free  field  and  the 
time  constant  in  seconds  (usually  0.4)  used  for  the  "inertia"  of 
the  receiver  in  dynamic  simulation.  The  program  will  then  gener- 
ate the  data  file  and  return  to  the  SWITCH  POINT. 

The  actual  ideal  surface  of  zero  COI  is  a hyperboloid  of  two 
sheets  whose  axis  of  rotation  is  parallel  to  the  z-axis  and  passes 
through  the  antenna.  For  comparison  purposes  it  is  convenient  to 
have  the  aircraft  travel  along  this  surface  and  see  how  the  real 
CDI  deviated  from  0.  If  the  glide  path  is  used  for  linear  flight 
in  the  near  field  (between  threshold  and  antenna)  large  CDI's  will 
occur  because  of  the  hyperbolic  shape.  The  program  will  allow  the 
user  to  generate  a hyperbola  which  is  the  intersection  of  the  0 
CDI  surface  and  the  plane  containing  the  runway  centerline  parallel 
to  the  z-axis.  To  do  this  the  user  types  a ’G’  for  the  flight  path 
type.  The  program  will  respond: 

INPUT  XO.  XF.  H: 

The  user  inputs  these  in  floating  point  free  field  format 
followed  by  a <CR>.  XO  is  the  x-coordinate  of  the  initial  receiver 
point  (Y0  is  zero  and  Z0  is  specified  by  the  hyperboloid),  XF  is 
the  x-coordinate  of  the  final  receiver  point.  H is  the  height  of 
the  main  carrier  element  in  the  antenna  array.  The  program  will 
then  respond: 

INPUT  * OF  POINTS,  VELOCITY,  TIME  CONSTANT 

These  are  input  as  above.  The  program  will  generate  the  file  and 
return  to  the  switch  point. 

4.4  SWITCH:  A 

This  switch  is  used  to  generate  the  antenna  description  file. 
The  program  will  respond: 
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INPUT  ANTENNA  FILE  NAME: 


The  user  inputs  the  5 character  file  name.  The  program  will  type: 
INPUT  ANTENNA  DESCRIPTION: 

The  user  inputs  a <40  character  antenna  description  to  be  used  as 
a plot  label.  The  program  will  then  type: 

INPUT  ELEMENT  VALUES: 


The  user  then  types  in,  in  free  field  floating  format,  a maximum 
of  8 fields  followed  by  a <CR>.  The  fields  have  the  following 
usage: 


field  » 
} 

2 

3 

4 
3 
6 

7 

8 


usage 

x-coordinate  of  clement  (usually  0) 
z-coordinate  r-£  element  (height) 
real  amplitude  of  carrier 
imaginary  amplitude  of  carrier 
real  amplitude  of  150  Hz  side  band 
imaginary  amplitude  of  150  Hz  side  band 
real  amplitude  of  90  Hz  side  band 
imaginary  plitude  of  90  Hz  side  band 


This  element  inputting  is  repeated  for  each  element.  After  the 
last  element  is  input  an  extra  carriage  return  is  typed.  No 
y- coordinate  is  input  for  the  elements.  This  is  because  nominally 
all  the  elements  have  the  same  y offset  (the  value  input  as  Y0 
under  switc.i:Y).  However,  a small  offset  correction  is  applied 
for  near  field  correction.  An  explanation  of  this  correction  may 
be  found  in  part  I (see  discussion  proceeding  Eq.  (33)).  This  is 
automatically  done  by  the  program.  As  the  first  clement  input  is 
assumed  to  have  the  correct  offset,  it  will  always  have  a value 
of  Y0.  Thus  the  main  carrier  element  should  be  input  first,  for 
example,  a null  reference  antenna  was  input  as  follows: 


a. 


INPUT  SWITCH:  Y 

INPUT  Y0,  LAMBDA:  300., 
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INPUT  SWITCH:  A 


5.  GLDPLT  PROGRAM  DESCRIPTION 


GLDPLT  is  a plotting  output  program  to  graph  the  CDI  informa- 
tion from  ILSGLD.  The  urer  runs  GLDPLT  which  then  types: 

INPUT  FILE  NAME  AND  AXIS  TYPE: 


The  user  then  types  in  a 5 character  (left- justified  and  blank 
filled  to  5 characters)  and  two  integer  fields  in  free  field 
format.  The  first  integer  is  the  switch  for  x-axis  type  and  the 
second  integer  is  the  switch  for  the  y-axis.  The  y-switch  has 
two  values,  1 and  2,  the  x-switch  has  three  1,  2,  and  3;  any  other 
values  will  terminate  the  program. 

The  switches  have  the  following  usage: 


y-switch*l 

y-switch*2 

x-switch*l 


x-switch-2 

x-switch*3 


this  plots  the  static  CDI  values 

this  plots  the  dynamic  CDI  values 

this  uses  the  altitude  angle,  in 
degrees  measured  from  the  origin  of 
the  receiver  point  as  x-coordinate 

this  uses  the  x-coordinate  of  the 
receiver  as  the  x-axis 

uses  the  time  in  seconds,  at  the 
receiver  point,  as  the  x-coordinate. 


After  the  input,  a plot  is  generated  and  the  program  returns  to 
the  input  and  asks  for  the  data  for  the  next  plot.  The  user  can 
give  the  same  file  name  to  plot  the  data  differently,  or  a new 
file  name  can  be  given.  This  would  be  done  when  multiple  runs 
were  done  before  plotting;  the  output  file  from  each  simulation 
run  carrying  its  own  name. 
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In  this  Appendix  a test  case  is  presented  which  the  programmer 
may  duplicate  to  gain  experience  in  running  the  model. 

The  test  case  uses  a Null  Reference  antenna  over  a terrain  as 
sketched  below. 


The  results  are  shown  in  Figure  9. 
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